Medium effects close to s- and p-wave Feshbach resonances in atomic Fermi gases 
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Many-body effects may influence properties, such as scattering parameters, nature of pairing, 
etc., close to a Feshbach resonance in the fermion BEC-BCS crossover problem. We study effects 
such as these using a tractable crossing-symmetric approach. This method allow us to include 
quantum fluctuations, such as, density, current, spin, spin-current and the higher-order fluctuations 
in a self-consistent fashion. The underlying fermion interaction is reflected in the "driving" term. 
We perform calculations here on both Bose-Einstein condensate (BEC) and BCS sides, and taking 
the driving term to be finite range, and of arbitrary strength. These are related to two-body singlet 
and triplet scattering parameters, and can be connected with experimental s— and p-wave Feshbach 
resonances. We include the £ = density and spin fluctuations, as well as I — 1 current and spin- 
current fluctuations. We calculate renormalized scattering amplitudes, pairing amplitudes, nature 
of pairing, etc., on both the BEC and BCS sides. We then compare our results qualitatively with 
experiments. 

PACS numbers: 71.10.Ay,34.50.-s, 05.30.Fk, 03.75.Ss 



I. INTRODUCTION 

The impact of ultracold atomic and molecular quan- 
tum gases on present-day physics is associated with the 
extraordinary degree of control that such systems offer to 
investigate the fundamental behavior of quantum mat- 
ter under various conditions Recent experimental 
achievements in the field of ultracold Fermi gases are 
based mainly on the possibility of tuning the scatter- 
ing length a s , in particular to values much large than 
the mean interatomic distance, by changing an external 
magnetic field [|J . This situation exists near the so-called 
Feshbach resonances. 

Resonances in general refer to the energy-dependent 
enhancement of interparticle scattering cross section due 
to the existence of a metastable state. For Feshbach res- 
onances, the metastable state is described in terms of 
coupling of a bound state of a subsystem to its environ- 
ment. By tuning a magnetic field, it is possible to obtain 
a quasidegeneracy between the relative energy of two col- 
liding atoms and that of a weakly bound molecular state. 
As the quasibound state passes through a threshold, the 
scattering length can be varied, in principle, from positive 
to negative infinity. The Feshbach resonances were ob- 
served in bosons [3H6[, in fermions between distinct spin 
states 043 j and in a single-component Fermi gas [loj . 
In this manner the interactions between the atoms can 
be strongly enhanced by an external magnetic bias field, 
giving rise to the BEC-BCS crossover phenomena [Til — 
[l3| . As a result of the atomic physics of the Feshbach 
resonance, the nature of the Cooper pairs in the BEC- 
BCS crossover is, however, not solely determined by the 
interaction strength or scattering length but, in principle, 
also depends on the width of the Feshbach resonance. In 
the limit of an infinitely broad resonance, the properties 



of the gases can be described by a single-channel theory 
that requires only the resonant scattering length as an 
experimental input. In general, however, a two-channel 
theory is needed. This is, in particular, true for the de- 
scription of the wave function of the Cooper pairs that 
plays an important role in the BEC-BCS crossover. 



Near a resonance, fluctuations can be quite important 
since the system is no longer in the dilute limit, and there- 
fore one may expect contributions from higher angular 
momenta and from quantum fluctuations or influence of 
collective modes. Generally, inclusion of many-body ef- 
fects in these systems is not easy, and one usually stops 
at the level of random phase approximation (RPA). 



In this paper, we aim at presenting a different scenario 
which emphasizes the role of the exchange particle-hole 
fluctuations in addition to RPA in the direct particle-hole 
channel. Our calculation is based on the induced inter- 
action model of Babu and Brown [3, [ll| , subsequently 
generalized and termed as the "crossing-symmetric ap- 
proach" . This takes into account a properly antisym- 
metrized, effective two-body interaction that reproduces 
the correct low-energy physics and also suppresses any 
spurious ground states. The physics is described in terms 
of the Landau interaction parameters and scattering am- 
plitudes to be discussed. In a previous study Gaudio et 
al. [2l| considered s-wave scattering by including I — 
density and spin fluctuation. We shall study s- and p- 
wave scattering and include both £ = density and spin 
fluctuations and £ = 1 current and spin-spin fluctuations. 



2 



II. THEORETICAL APPROACH 
A. The crossing-symmetric method 

The crossing-symmetric method [I3 - [l6j was formu- 
lated to calculate the effective quasiparticle interactions 
in Fermi systems. Due to an appropriate compromise 
between microscopic and phenomenological approaches, 
it has been successfully applied to a number of Fermi 



systems: liquid 1 He 



any applied 

EMI [13, 



paramagnetic metals 
nuclear matter [13 
gases pH 



'He- 4 He mixtures 



161 ]. heavy- fermion systems 
and ultracold atomic Fermi 
It has been known [II [H [H, [H that a 



consistent Fermi-liquid theory cannot be formulated in 
terms of short-range effective interactions alone; collec- 
tive excitations generated by these must be exchanged 
between quasiparticles. 



the Lindhard functions or density-density and current- 
current correlation functions, respectively. The first term 
in Eqs. (|3]) and Q is the so-called direct interaction. The 
direct term is designed to convey the fact that two quasi- 
particles can directly scatter via some effective potential 
and repeatedly so , as in a T-matrix. It is of short range 
and contains information about the underlying Hamilto- 
nian of the system under consideration. Thus, it is the 
"driving" term. The induced term is of somewhat longer 
range since two particles can scatter via an interaction 
mediated by another particle. Equations. (O and (d]) are 
nonlinear coupled equations. To solve these, we need to 
do Legendre projections: 



K£ = ^2F l s > a p l (cos6 L ), 



pp' 



pp 



Df' a Pi(co80 L ) 



(5) 
(6) 



The main point is that the contributions to Landau 
interaction /£5,' can be separated into two parts [l5j |: 



aa 
PP' 



«pp' 



PP' 



J P 



(1) 



where the induced part I™, , a function of the Landau 
interactions f°Z, themselves, is particle-hole reducible 
in the exchange particle-hole (u) channel, whereas 
the direct part d^, is not particle-hole reducible in 
either the direct particle-hole (i) channel or the crossed 
particle-hole (u) channel. It is important to note that 
the direct interaction is model dependent, as it gives 
information about the underlying Hamiltonian, and 
that the induced interaction is a purely quantum effect, 
arising from the exchange diagrams required to antisym- 
metrize the effective two-body scattering amplitude. 

For a two-component fermionic system, the Landau 
interaction can be expressed as 



Ainsworth et al. [16j treated Dq, Dq, and Df phenomeno- 
logically so as to reproduce the empirical Landau param- 



eters F£, Fq, F{, and F" and predicted the higher-order 



obtained 




Then effective pairing interaction can be 



f f = [A s - 3A a - (A{ - 3AI)] /4 = A s /4, 
[AZ+A$-(Al +A?)]/12 = A t /12, 



eff 
9t 



where 



A s.a 



Ff 



1 + F Z S '7(2/ + 1) 



(7) 
(8) 

(9) 



is the Z-partial wave symmetric (s), antisymmetric (a) 
scattering amplitude and A s and A t are the pairing- 
channel singlet and triplet scattering amplitude respec- 
tively. The vanishing of forward scattering of two parti- 
cles of equal spin yields the Landau sum rule ^2i(Af + 
Af) = 0, which provides a test for obtained Landau pa- 
rameters. 



F aa , = F s , 
pp PP 



(2) 



B. The driving term near a Feshbach resonance 



where F, 



pp' 



and 



pa 
pp' 



can be related to induced- 



interaction equations as follows [If 

1 F°U (q')F£ 
21 + FgUoffl 
FfUx^Ff 



F s , = D s , 
pp' pp' 



1 



44 



3 FgU (q')Fg 
2 1 + F«U (q>) 
3 F?U x {q>)F? 



= D a , 
pp 



1 - 



4k 2 F 



1 + FfU^q') l + FfUttf) 

1 fggoCgQfg 1 FgU {q')FS 
21 + F$U (q>) 
FtU^q'Wl 



(3) 



2 1 + F-Uo(q') 



(4) 



1 + FfUtiq') l + F?U 1 (q>)_ 
The momentum transfer in the crossed particle-hole 



channel is q' = 
dau angle 9l 



cos 



= kp\J\ — cos 0^, with the Lan- 
(p-p'). Uo(q') and Ux{q') are 



In the crossing-symmetric approach, the form of the 
direct interaction used to derive the induced interactions 
must be determined. In general, it represents the sum of 
all particle hole irreducible interactions. A self-consistent 
calculation could be performed starting with a T-matrix 
direct interaction if a more general interaction were used 
to derive the T-matrix calculation. According to the pro- 
posal by Bedell and Ainsworth [25[ , the direct interaction 
is the Fourier transform of an effective quasi-particle po- 
tential. From this potential the quasiparticle scattering 
amplitude /k(0) is given by 



H4>) 



4tt 



e^ r V eff (r,k)d 3 r, 



(10) 



where ft = 1, q 2 = |k-k'| 2 = 2fc 2 (l-cos (f>), 2k 2 = k 2 F (l- 
cos#) and the quasiparticle mass m* — m(l+F-f/3). The 
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relative momentum of the incoming (outgoing) particles 
is k(k'), and the angle between the incoming and scat- 
tering plane is <f>. Equation (fT0|) is restricted to the Fermi 
surface; therefore /k(</>) depends on only two variables, 9 
and 4>. According to the effective range expansion, the 
effective potential 14// (r, k) is, in general, nonlocal and 
can be expanded in powers of k 2 , 

V eff (r, k) = U(r) + ±k 2 r 2 W(r) + (11) 

where U(r) and W(r) are local potentials. Keeping or- 
der of k to 2 yields a three-parameter approximation to 

AW, 

f k ((j>) ~ — [-a, + 3k 2 a t (l - cos<j}) - 3k 2 & t l , (12) 
to 

where a s = m f r 2 U(r)d 3 r, at = (m/9) J r 4 [/(r)d 3 r, and 
b t = (m/9) Jr 4 W(r)d 3 r. As pointed out in Ref. [H, the 
triplet quasiparticle scattering volume at and the nonlo- 
cal part of the effective potential bt are p wave in nature 
and thus sample relatively little of the repulsive core of 
the bare interaction. Therefore they should not have a 
strong density dependence. The a t is a finite-range cor- 
rection to the contact interaction; this allows some inter- 
action between particles of the same spin. The nonlocal 
piece of the direct interaction results from the coupling 
of quasiparticle currents. The direct interaction for par- 
ticles of parallel spin is 

Air 

d^{e,<t>) = — -\M<i>)-M<t> + -«)] 

to* 
12-k 

= kpaAl — cos 9) cos (j). (13) 

to 

Similarly for particles of opposite spin, 

4-7T 

^(0,0) = ~M<f>) 

to* 

— — \a s — —kpdtil — cos#)(l — cos0) 
m 2 

+ ^4&t(l-cos0)]. (14) 

The standard Landau parameters are the q — values 
of momentum dependent functions i^ s ' a (q). In the limit 
of q — > 0, cos9 = cos6*i, cos</> = 1. dP~, d^- 1 - could be 
rewritten as 

n ' 4 f ^(P • P') = —k 2 F a t (l - cos0 L ), 



pp TO 

dp P > = 4 l + 4 l Pi(P ■ P') = ^[a. + \k 2 F b t {\ - cos9 L )] 



Solving the above equations, we obtain 



d n 



12tt 

TO 

4tt 



7 2 

KpOf> 



+ —kibt, 



12tt, 



m 

TO 



-k F a t , 
kpb t . 



(15) 



The direct interactions D s and D a are linear combina- 
tions of (ftt and d^: 

DS = NM {d n + d n^ (16) 



(17) 



where N(0) = kpm* /it 2 is the density of states at Fermi 
surface. For a pure s-wave resonance, we set a t = 0, bt = 
0: 

D s = —a s N(Q) = U /2, D{ = (18) 

TO 

= - — a s N(0) = -U /2, D" = 0. (19) 

TO 

For a pure p-wave resonance, we set a s = 0: 
lirk 2 

D s = ^(2a t + bt)N(0),Df = -D°; (20) 

TO 

D- = ^(2a t - 6 t )iV(0), Dt = -D%. (21) 



C. Connection to the usual scattering parameters 

The two-body scattering amplitude is 



i=0 



J2(2l + l)fi(k)Pi(cos( 



(22) 



2=0 



where the Ith angular momentum channel is given by [2 

A: 2 ' 



fi(k) 



+ r t k 2 - ik 2l+1 



(23) 



According to the effective range expansion, we have 



1 ,2 

cot rf i = h rik + 

ai 



(24) 



From the induced interaction description, as q = 0, 2fc 

fp ( 



kp(l — cos 9). Equation (TT2")) becomes 



to* 3 3 

fk(P) = —[-as ~ ^kpb t + -k 2 F b t cos6] (25) 
to 2 2 

Combining Eqs. (j2"2"l) - (|2"5"|) . we reach 
k 2 



TO* 1 

~^2 



k 2 F b t 



-dj 1 + r\k 2 — ik 3 



yj{k 3 cot ?7i) 2 + k 6 ' 
(26) 
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FIG. 1. (Color online) The s-wave scattering length as a 
function of the magnetic field B on both sides of the Feshbach 
resonance in 40 zK, using data from Regal and Jin [27|]. The 
density here is n = 5.8 x 10 13 cm~ 3 , with AB = 9.7G, and the 
Feshbach resonance occurs at Bo = 224. 21G. The scattering 
length is measured in terms of the Bohr radius. 
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FIG. 2. (Color online) The scattering cross section of a s- 
wave resonance between 40 K atoms in \f = 9/2, mj = —9/2 > 
and |/ = 9/2, m/ = —7/2 > states. The data used [10J for 
calculation are as follows: number density is n p k = 1.5 x 
10 13 cm -3 , with AB = 7.8G, and the Feshbach resonance 
occurs at B = 202. 1G. 



III. CALCULATIONS AND RESULTS 

A. Close to an s-wave resonance 

In the vicinity of a Feshbach resonance, the s-wave 
scattering length a(B) is described approximately by [28[ 



a(B) = a bg 1 - 



AB 

B-Br 



(27) 



where B is the applied field, a{, g is the background scat- 
tering length, and Bq is the field at which the resonance 
occurs. The resonance width AB is proportional to the 
strength of the coupling between the open and close chan- 
nels. For 40 K, a^g — 174ao [29], where ao is the Bohr 
radius. Based on this, we construct the driving terms 
and solve Equations ^ and We find that the scat- 
tering length tends to smooth out as it approaches the 
resonance, as shown in Fig. [TJ On the BCS side far from 
the resonance, the many-body results give exactly two- 
body physics. When the system is driven to the reso- 
nance, two-body scattering gives a diverging unrenormal- 
ized (bare) scattering length. Many-body exchange fluc- 
tuations greatly suppress the divergence of the scattering 
length. In this region, the exchange fluctuations are quite 
strong and act as a feedback to the system. On the Bose- 
Einstein condensate (BEC) side far from the resonance, 
the medium effects reduce the effect of the interaction. 
This is similar to what occurs on the BCS side far from 
the resonance. The suppression of divergence close to the 
resonance is also suggested in the many-body renormal- 
ized scattering cross section, as shown in Fig. [5J where 
scattering occurs between two hyperfine species, namely 
|/ = 9/2, to/ = -9/2 > and \f = 9/2, m/ = -7/2 >. 
Here / is the total angular momentum, and to/ is the 
corresponding magnetic quantum number. As can be 
seen, away from the resonance, the medium effects are 
small. However, close to resonance, the medium effects 
strongly modify the scattering cross section. 



B. Close to a p-wave resonance 

A p-wave resonance is distinct from an s-wave (I = 0) 
resonance in that the atoms must overcome a centrifugal 
barrier to couple to the bound state. It is sensitive to 
the temperature and the magnetic field. We need two 
parameters to characterize the p-wave resonance: scat- 
tering volume at and effective range n. The magnetic 
field dependence of a t and r\ is obtained from the fit- 
ting formula given by Ticknor et al. [30]. Therefore we 
can construct the driving terms and solve equations ([3]) 
and (UJ). The p-wave resonance could be tuned by the 
magnetic field to occur between two atoms in the hy- 
perfine states of ]/, to,/ >= 1 9/2, — 9/2) and ]/, to./ >= 
1 9/2 , — 7/2). The joint state of the atom pair will be 
written as 1/iTO/J |/2TO/ 2 ) |^to^). The many-body effects 
in a p wave are less severe than in an s wave, as shown 
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FIG. 3. (Color online) The scattering cross section of a p- 
wave resonance between 40 K atom pairs in the state \f = 
9/2, m/ = -9/2)|/ = 9/2,771/ = -7/2 > \£ = l,m e = 0) at 
T = 3.2[iK. 
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FIG. 4. (Color online) The many-body scattering cross sec- 
tion of a p-wave resonance for 40 K atom pairs in the |/ = 
9/2,771/ = -9/2)1/ = 9/2,m/ = -7/2) \£ = l,m t = 0) state 
at various temperatures. 
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FIG. 5. (Color online) The singlet and triplet scattering am- 
plitudes (A s and At) for a p-wave resonance between 40 K 
atom pairs in the state \f = 9/2, m/ = — 9/2) |/ = 9/2, m/ = 
-7/2}\£ = l,mt = 0) at T = O.lfxK. 



T = 5.0/iK (green dots) shows a double-peak feature re- 
sulting from the strong energy dependence of the cross 
section. The pairing-channel scattering amplitudes (A s 
and At) are shown in Fig. [5j The triplet scattering ampli- 
tude A t is negative on the BCS side, indicating a p-wave 
supcrfluid pairing. On the BEC side, the singlet scat- 
tering amplitude is negative, indicating the formation of 
s-wave molecules. The p-wave Feshbach resonances of- 
fer a means to experimentally study anisotropic inter- 
actions in systems other than identical fermions. On 
resonance the p-wave cross section becomes comparable 
to the background s-wave scattering. This means that 
it could have an equally important role in determining 
the collisional behavior and mean-field interaction of the 
quantum gases. Utilizing p-wave Feshbach resonances, 
the Joint Institute for Laboratory Astrophysics (JILA) 
group has successfully created p-wave molecules [31( . 



in Fig. [3] In Fig. [31 we plot the scattering cross section 
for a p-wave resonance for 40 K atom pairs in the state 
|/ = 9/2, m/ = —7/2) |^ = l,m e = 0). Here t = 1 
is the orbital angular momentum quantum number, and 
Tfi£ — is the corresponding magnetic quantum num- 
ber. Away from the resonance, the many-body correc- 
tions due to medium effects are small. Distinct correc- 
tion only appears at a region close to the Feshbach res- 
onance. Figure 0] shows the temperature dependence of 
the many-body scattering cross section. We can see that 
the lower the temperature is, the higher the resonance 
peak is. As the temperature rises, the resonance cross 
section broadens. The position of the resonance changes 
slightly with temperature due to the temperature depen- 
dence of the scattering cross section. The resonance at 



C. Close to a p-wave resonance with an s-wave 
background 

For a p-wave resonance, the JILA group found that 
there exists non-negligible off-resonant scattering in the 
ultracold Fermi gas of 40 K. Our model can be employed 
to study the situation when the singlet correlation (s- 
wave scattering) and the triplet correlation (p-wave scat- 
tering) are both present. This situation can occur when 
the resonant magnetic field for p waves and s waves are 
well separated. By tuning the magnetic field around p- 
wave resonance one can realize a p-wave resonance with 
an s-wave background. The driving terms can then be 



6 




200 




FIG. 6. (Color online) The singlet scattering amplitude A 3 
for a p-wave resonance: (a) without s-wave background and 
(b) with s-wave background where Uo = —1.5. 
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FIG. 7. (Color online) The triplet scattering amplitude At 
for a p-wave resonance: (a) without s-wave background; (b) 
with s-wave background where Uo = —1.5. 



modeled as 



3tt. 



2tt 



— kp(2a t + b t )N(0) + —a s N(0), 
m m 



3tt. 



2tt 



Do = — 4(2o* - b t )N(0) a,N(0), 

m m 

37T 

D{ = 4(2ot + bt)N(0), 

m 

D? = k 2 F (2a t -b t )N(Q). 

m 



(28) 



Here s-wave background interaction is characterized by 
Uo = ^-a s N(0). By fixing U and varying the parame- 
ters l/kpa t and 1/kpbt, we can evaluate the pairing chan- 
nel scattering amplitudes (A s and A t ) in full parameters 
space through solving crossing-symmetric equations ([3]) 
and ([4]). From the obtained Landau parameters Ff' a , by 
s-p approximation (24)], one can straightforwardly con- 
struct pairing-channel scattering amplitudes via 



Asinglet — Aq — 3j4q — ^1 "F 3^4 1 , 
^triplet = "F ^0 — ^1 — ^l" 



(29) 



The calculated singlet scattering amplitude is plotted in 
Fig. [S] When there is no background s-wave scattering, 
the singlet scattering amplitude becomes very large when 
it crosses the unitary limit (strongly interacting regime) 
of either parameter l/k F a t or l/k F b t , as suggested in 
Fig. \§la). To investigate the effects of the s-wave back- 
ground on the p-wave resonance, we choose Uq = —1.5 
for illustration such that it has a noticeable effect. In 
the presence of background s-wave scattering, the singlet 
scattering amplitude becomes negative and shows a pro- 
nounced peak when it crosses the unitary limit of both 
parameters 1/kpat and 1/kpbt- Away from the unitar- 
ity limit, the singlet scattering amplitude A s is mainly 
negative, manifesting trends of background s-wave pair- 
ing. The purpose of introducing Uq is to introduce back- 
ground singlet pairing interaction. The competition be- 
tween the p-wave scattering and background s-wave scat- 
tering determines the sign of the singlet scattering ampli- 
tude. When the p-wave interaction is weak, the proper- 
ties of the system are dominated by s-wave behavior with 
a constant negative singlet pairing amplitude, as can be 
seen in the Fig. \6lb). 

To investigate p-wave pairing, we plot the triplet scat- 
tering amplitude A t in Fig. [JJ Interestingly, the triplet 
scattering amplitude A t is sensitive to parameter 1/kpat 
but insensitive to parameter 1/kpbt except building up 
a bump near the unitarity regime in the presence of 
background s-wave scattering. Remarkably, in the ab- 
sence of background s-wave scattering A s is antisymmet- 
ric with respect to parameter 1/kpat while symmetric 
with respect to parameter 1/kpbt- These properties are 
closely related to the parametrization and symmetrical 
features in Equation (28), where s-wave scattering only 
contributes to the £ — channel in the driving terms. 
The background s-wave scattering shows its existence by 
shifting the triplet scattering amplitude in total toward 
the negative side. In addition, it enhances the effect of 
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parameter 1/kpbt, especially near resonance. The anti- 
symmetry of A t with respect to parameter \/k\.a t is lost; 
in contrast, the symmetry of A t with respect to param- 
eter l/kpb t is still preserved. It is necessary to point 
out that when both A s and A t are negative, the actual 
pairing symmetry depends on their relative magnitude. 



IV. CONCLUDING REMARKS 

By using the crossing-symmetric method to treat 
many-body Fermi systems, one generally solves nonlin- 
ear coupled crossing-symmetric equations for four-point 
vertex functions in particle-particle, particle-hole, and 
exchange particle-hole channels. In appropriate limits 
on the Fermi surfaces, these vertex functions become 
the Landau quasiparticle interaction F(q) and the 
scattering amplitude A(q). For isotropic systems, these 
can be expressed in Legendre polynomials giving Landau 
parameters Ff and F" in spin-symmetric (s) and spin- 
antisymmetric (a) channels. From the obtained Landau 
parameters, one can calculate several thermodynamic, 
transport, and pairing properties of a system. 



We study s- and p-wave Feshbach resonance with the 
crossing-symmetric method. Our findings are as follows: 
(1) many-body exchange effects may be important close 
to a Feshbach resonance. Renormalized physical quan- 
tities get smoothed out at the Feshbach resonance. In 
particular, we find that the particle-hole exchange fluc- 
tuations introduce an effective scattering length which 
has been substantially reduced close to resonances. (2) 
For a p-wave resonance, the triplet scattering amplitude 
is negative on the BCS side, indicating a triplet (p- wave) 
supcrfluid pairing. (3) Background off-resonant scatter- 
ing has some effects on the singlet and triplet scattering 
amplitudes, which may influence the pairing symmetry 
of the ground state. 

ACKNOWLEDGMENTS 

We are grateful for stimulating discussions with K. S. 
Bedell and S. Gaudio. This work was supported by the 
NSFC under Grants No. 10934010 and No. 60978019, 
the NKBRSFC under Grants No. 2009CB930701, 
No. 2010CB922904, No. 2011CB921502, and No. 
2012CB821300, and the NSFC-RGC under Grants No. 
11061160490 and No. 1386-N-HKU748/10. 



C. Chin, R. Grimm, P. Julienne, and E. Tiesinga, Rev. 
Mod. Phys. 82, 1225 (2010) [16 
S. Giorgini, L. P. Pitaevskii, and S. Stringari, Rev. Mod. 
Phys. 80, 1216 (2008) [17 
S. Inouye, M. Andrews, J. Stenger, H. Miesner, D. M. 
Stamper, and W. Ketterle, Nature 392, 151 (1998) [18 
P. Courteille, R. S. Freeland, D. J. Heinzen, F. A. van 
Abeelen, and B. J. Verhaar, Phys. Rev. Lett. 81, 69 [19 
(1998) 

J. L. Roberts, N. R. Claussen, J. P. Burke, C. H. Greene, [20 
E. A. Cornell, and C. E. Wieman, Phys. Rev. Lett. 81, [21 
5109 (1998) 

A. Marte, T. Volz, J. Schuster, S. Durr, G. Rempe, [22 
E. G. M. van Kempen, and B. J. Verhaar, Phys. Rev. 
Lett. 89, 283202 (2002) [23 
T. Loftus, C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. 
Jin, Phys. Rev. Lett. 88, 173201 (2002) [24 
K. M. O'Hara, S. L. Hemmer, S. R. Granade, M. E. 
Gehm, J. E. Thomas, V. Venturi, E. Tiesinga, and C. J. [25 
Williams, Phys. Rev. A 66, 041401 (2002) 

K. Dieckmann, C. A. Stan, S. Gupta, Z. Hadzibabic, [26 
C. H. Schunck, and W. Ketterle, Phys. Rev. Lett. 89, 
203201 (2002) [27 

C. A. Regal, C. Ticknor, J. L. Bohn, and D. S. Jin, Phys. 
Rev. Lett. 90, 053201 (2003) [28 

D. M. Eagles, Phys. Rev. 186, 456 (1969) 

A.J. Leggett, Modern Trends in the Theory of Condensed 

Matter (Springer- Verlag, Berlin, 1980) [29 

P. Nozieres and S. Schmitt-Rink, J. of Low Temp. Phys. 

59, 195 (1985) ' [30 

S. Babu and G. E. Brown, Ann. Phys. 78, 1 (1973) 

K. F. Quader, K. S. Bedell, and G. E. Brown, Phys. Rev. [31 



B 36, 156 (1987) 

T. L. Ainsworth, K. S. Bedell, G. E. Brown, and K. F. 
Quader, J. of Low Temp. Phys. 50, 317 (1983) 
K. F. Quader and K. S. Bedell, J. of Low Temp. Phys. 
58, 89 (1985) 

L. Yi, K. S. Bedell, and T. L. Ainsworth, Bull. Am. Phys. 
Soc. 28, 655 (1983) 

K. S. Bedell and K. F. Quader, Phys. Rev. B 32, 3296 
(1985) 

O. Sjoberg, Ann. Phys. 78, 39 (1973) 

S. Gaudio, J. Jackiewicz, and K. S. Bedell, Phil. Mag. 

Letts. 87, 713 (2007) 

K. S. Bedell and K. F. Quader, Phys. Lett. 96A, 91 
(1983) 

K. S. Bedell and K. F. Quader, Phys. Rev. B 30, 2894 
(1984) 

B. R. Patton and Z. Zaringhalam, Physics Letters 55A, 
95 (1975) 

T. L. Ainsworth and K. S. Bedell, Phys. Rev. B 35, 8425 
(1987) 

L. D. Landau and E. M. Lifshitz, Quantum Mechanics 
(Permagon, Oxford, 1994) 

C. A. Regal and D. S. Jin, Phys. Rev. Lett. 90, 230404 
(2003) 

S. J. J. M. F. Kokkelmans, J. N. Milstein, M. L. Chiofalo, 
R. Walser, and M. J. Holland, Phys. Rev. A 65, 053617 
(2002) 

M. Greiner, C. A. Regal, and D. S. Jin, Nature 426, 537 
(2003) 

C. Ticknor, C. A. Regal, D. S. Jin, and J. L. Bohn, Phys. 
Rev. A 69, 042712 (2004) 

J. P. Gaebler, J. T. Stewart, J. L. Bohn, and D. S. Jin, 



Phys. Rev. Lett. 98, 200403 (2007) 



